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Abstract 

Recently, a functional integral representation was proposed by Weller , in which 
the fermionic fields strictly satisfy the constraint of no double occupancy at each lattice 
site. This is achieved by introducing spin dependent Bose fields. The functional integral 
method is applied to the single impurity Anderson model both in the Kondo and mixed- 
valence regime. The f-electron Green's function and susceptibility are calculated using 
an Ising-like representation for the Bose fields. We discuss the difficulty to extract 
a spectral function from the knowledge of the imaginary time Green's function. The 
results are compared with NCA calculations. 

keywords: functional integral, single impurity Anderson model, f-electron Green's 
function 



1 Introduction 



In systems of correlated fermions on a lattice one often encounters the presence of a very 
high local Coulomb repulsion U between two particles on the same lattice site. In this case 
simple perturbation theory in the parameter U no longer provides a good approximation. In 
the limit U — oo, the interaction term in the Hamiltonian can be eliminated completely if a 
constraint is introduced which allows only empty and singly occupied states. 

Functional integral techniques are a very powerful tool for the investigation of lattice fermion 
systems and there are several ways to incorporate the constraint in the functional integral. 
Starting for example with a slave boson approach, the constraint is guaranteed by delta- 
functions in the integration measure. 

The method we consider here is based on the idea of projecting out the doubly occupied 
sites already in the derivation of a coherent state functional integral. The new theory was 
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proposed and studied by Weller |I| and has already been applied to the one-j| and two-|3|. Q 
dimensional Hubbard model. 

In this paper, we focus on the single impurity Anderson model (SIAM) with infinite local 
Coulomb repulsion at the impurity site. This model was proposed by Anderson || to describe 
the properties of metals containing magnetic impurities and has been intensively studied (see 
e.g. |], [7|, P|), so that good approximate results are available today. However, these methods 
do not work for all parameter ranges of the model or cannot easily be extended to the 
lattice case. Therefore there is still some need for alternative methods to solve the impurity 
Anderson model. 

We have to state here that also with our method, we do not succeed in obtaining more 
satisfying results for the SIAM. Nevertheless, we think it is worth publishing the results 
because the new method is conceptionally quite different to the approaches cited above and 
on the other hand may be useful in the future if the limitations by computer time become less 
serious. Even analytical approaches based on this functional integral method may become 
possible, perhaps also in the investigation of other models of strongly correlated electron 
systems. 

The rest of the paper is organized as follows. After an introduction to the new functional 
integral method (section 2) we derive an effective action for the impurity states of the SIAM 
by integrating over the conduction electron degrees of freedom (section 3). This action is 
used as the starting point for the numerical investigation described in section 4. In section 5 
we present results for the imaginary time Green's function G(t) from which the f-occupancy 
can be derived immediately. Using a fit procedure we obtain the spectral function A(u). To 
demonstrate the limitation of the fit procedure we show that various quite different spectral 
functions can result in very excellent fits to the same given G(t). We also discuss the case of 
the impurity f-level lying at the Fermi level. Here we expect the spectral function to consist 
of a single Lorentzian. Within this restriction we get a broadening and renormalization 
of the f-level consistent with calculations within a second order perturbation theory in the 
hybridization. We compare the results with NCA (Non Crossing Approximation) calculations 
for both parameter sets. With an additional magnetic field coupling to the impurity state 
we calculate the static susceptibility and the effective magnetic moment. 



2 Theory 

In this section we follow closely the method originally published by Weller JE]. We want 
to introduce the functional integral method by considering first the simple example of a 
fermionic system with one site and infinite Coulomb repulsion 

H^etflU + Uflfiflft (U^oc) . (1) 
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Introducing Hubbard operators 



X 0a = \0)(a\ and X aO = \a)(0\ (2) 
which project out the doubly occupied states, the Hamiltonian can then be written as 

H = }^ e{X a0 X 0cr . (3) 



Essential for the calculation of statistical properties of the fermionic system at finite tem- 
peratures is the partition function 



Tr 



knT' 



(4) 



with the trace in the restricted Hilbert space. We divide exp(—/3H) into a product of N equal 
operators exp(— AtH), thereby defining steps r n = (n — l)Ar (n = 1, . . . N, (At = (3/N)) 
in the imaginary time interval [0,(3]. At each time step we insert unity operators l n which 
project out the doubly occupied states. 



Tr 
Tr 



e -ArH e -ArH _ _ _ ^AtH 



-AtH-h ^-ArH-tr 

JIat-1 



l 2 e 



-AtH 



with 



o)(o| + |T)(T| + ||)U 



(5) 

(6) 
(7) 



The unity operators actually do not depend on n. This index is used to distinguish auxiliary 
fields which will be introduced at each time step n (see (8)). The most elegant way to derive 
a functional integral for fermionic systems is to use projectors expressed in terms of coherent 
states (see e.g. M). Following Weller M, we make the ansatz: 



I, 



dipld^ n \ip n ,b n )(ip n ,b n 



with the coherent states 



0) = exp 



which are constructed to be eigenstates of the operator X { 



Oct 



X 0rT | ip n , b n ) = %l) n b nu I if> n , b n ) 



(8) 



0) (9) 



(10) 



The ip are Grassmann variables obeying the usual anticommutation relations. The 6- variables 
are ordinary complex numbers and do not have a direct physical meaning. They just count 
the spin multiplicity for each step in the imaginary time interval. Because of the restriction 
to zero or singly occupied sites and because we used spin carrying 6-variables in the ansatz 
(9), we only need one Grassmann field for each time step n. 
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From requiring the equality of equations (6) and (7), it is straightforward to prove that the 
following equations have to be fulfilled: 



= 1 ; 



b b ?1 = o; yi = i 



'11} 



lb 1 Jb * Jb 

In JE] several possible representations for the 6-variables are presented. Below we will use a 
special spin representation (see (30)). 

Replacing the trace by 



Tr\ 



d^ld^i(^i,&i| . . . | ipi, h) 



(12) 



the next step is to evaluate the matrix elements (i/; n , b n \e~ ArH \ ip n -i, b n -i). For nontrivial 
Hamiltonians this can only be performed in the limit At — > 0, and therefore the calculations 
become exact only in the limit iV — > oo. 



Finally we obtain the functional integral representation for the partition function: 

Z = JJ V^Vi) exp i)\ b) 

with: 



(13) 



N 



Vrn-1 



(14) 



n=l 



and 



Zty+Zty = y d^jvd^TvdV'jv-idV'iv-i • • • d^Id^i (15) 

Note that — due to the J2a b* av b an -\ term — we do not arrive at the usual expression with 
the time derivative (ip^-^ip) in the kinetic part of the action! 

The functional integral representation for the imaginary time Green's function 



- G(t) = (X a0 (r)X 0a (0)} = ^Tr [e^ H X a0 (r)X 0a (0) 



(16) 



can be derived in a quite similar way as described for the partition function. The result is 
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(17) 



with 



T = (n-1) — , n = l...N 

v N 



Note that in the integrand Grassmann variables always appear in products with 6-variables. 
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3 Action for the (U = oo) Single Impurity Anderson 
Model 



Now we want to derive an effective action for the f-electrons of the SIAM by integrating out 
the conduction electron degrees of freedom. This action will be the starting point for the 
numerical investigation described in section 4. 

In the limit U — > oo the Hamiltonian of the SIAM is given by: 



a ka 

+ E v ( x *oCk* + c\ a X ( 



ka 



(19) 



The operators c\ a (c^) create (annihilate) a conduction electron with spin a and wavevector 
k. The last term describes the hybridization between the localized f-state and the conduction 
electrons with the hybridization matrix element V independent of k. The energies Ek and £f 
are to be measured from the Fermi level. 



The action corresponding to the above Hamiltonian is: 



N ( 8 

n=l V iV a 

+ Y.Xkanii 1 

nka \ 



A. 



I 

N 



£k)~X.kan—l X,kan 



N 



an—l 



nka 



(20) 



For the conduction electron operators c\. a and Cfc CT (whose corresponding Hilbert space is not 
restricted) we applied the standard method of replacing operators by Grassmann variables 
(see e.g. 0). Now we have to integrate over three types of fields: the complex fe-fields, the 
Grassmann variables ip for the f-electrons and the Grassmann variables x f° r the conduction 
electrons. The partition function for the whole system then reads: 

V^Vi> J V X ] V X exp [S(ij } V f , X, x\ b)\ (21) 

The integrations can in principle be performed in several ways; we choose to first integrate 
over the Grassmann fields of the conduction electrons using the formula: 



/ v x ] v x 



cxp 



(det H) exp 



Eei("~%< 



(22) 



E x\H ijXj + Xitl + E xUi 

ij i i 

The resulting action 

S eff (^, V + , 6) = Sfty, ^, b) + (^-) Yl iiC+ii-iL-i Y.( H ~ l )nm{k) (23) 

\ / anm k 
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(Si is the first term at the right hand side of (20).) reduces the problem to that of a sin- 
gle electron with a time dependent coupling mediated by the conduction electron Green's 
function if -1 with the matrix elements: 



(H )nm(k) 



Using the formula 



(l + e _/3efc ) X exp(( 



n> m 



1 + e ^ efc ) exp ( (N — (m — n))jrEk) '■ n<m 



with 



Uip^Uip exp 







detM 



M nm (b) = (1 - j^ef)(£bl n b cm -. 1 )8n-i,m - S nm 



+ 



'/3V 



} a k 

the integration over the remaining Grassmann variables V'l.jv leads to 

Z = (det H) f det M(b) 

Jb 

for the partition function. For the imaginary time Green's function we get 

- G(t) = (deti?)i / 6 CTl 6; n (M- 1 (6)) 1 , n detM(6) . 
Z Jb 



(24) 



(25) 



(26) 



(27) 



(28) 



Note that this equation defines G(r) only for the r-values r = (n — (n = 1, . . . N) . 
G(P) cannot be calculated directly within this approach but is simply related to G(0) by 

G(fi) = 1 - 2G(0) (29) 

in the case of no double occupancy. We do not actually need to calculate det H because this 
factor cancels in the Green's function. 



4 Numerical Investigation 

The resemblance of the theory to slave boson techniques and the success of slave boson mean 
field approaches may suggest that also in our case a saddle-point like approximation for the 
6-fields is possible. But due to the unusual kinetic part in the action any replacement of 
e -g- J2a b* an b a n-i by a constant number not equal to one leads to divergences in the limit 
N —>■ oo. So far we did not find any analytical approach with the action (23) as a starting 
point and therefore restrict ourselves to numerical results in this paper. 

The most suitable representation of the 6-fields for a numerical implementation of the equa- 
tions (27) and (28) is an Ising-like representation 

hn = 1 h n = ±1 (30) 
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/••• (31) 

so that the sums in (27) and (28) contain 2 N terms. For each of these contributions we have 
to calculate the Matrix M defined by (26). The sum over k in (26) is calculated using a 
constant density of states for the conduction electrons in the interval [— D,D\. 

Y,{H- l ) nm {k) = J depie)^- 1 )^) = p f D D &E{H- l ) nm {e) (32) 

These integrations are independent of the b fields and therefore have to be performed only 
once. NAG routines are used to calculate the inverse and the determinant of M. For N = 20 
the computation and summation of all contributions to G(t) takes about one day CPU time 
on a workstation. 

At this point, one would like to reduce the computer time by restricting the summation to 
the most important contributions. This idea fails in our case! Due to the J2a ^Cn^n-i-term, 
contributions to the Green's functions can differ by a factor 2° to 2^, so that simple Monte 
Carlo methods (like the Metropolis algorithm) are unable to explore the whole phase space 
of the 6-variables. On the other hand we experienced, that taking into account only the most 
important contributions, the result is far away from the exact result received by summing 
over all contributions. 



5 Results 

Fig. 1 shows the imaginary time Green's function for the parameters E{ = — 0.2,V = 0.22 
and the inverse temperature (3 = 3 (energies in arbitrary units). It also contains the result 
of an NCA calculation which will be discussed below. 

The f-occupancy equals 

m = n T + ni = 2 ■ G(t = 0) = 0.72 . (33) 

Although G(t) seems to be rather structureless and indeed has no direct physical meaning, 
it is related to the spectral function by the transformation 

_ G(t) = r due-" \ A<yU \ . (34) 

It is easy to calculate G(t) with the knowlegde of A(u) but not vice versa. A lot of methods 
can be found in the literature to overcome this problem (see e.g. |IJ) but determination of 



the spectral function A(u) out of G(t) remains an extremely ill posed problem as will be 
shown below in an example. 

In order to allow the occurrence of an Abrikosov-Suhl resonance in addition to a broad peak 
at the f-level position we tried to fit the numerical G(t) data with a superposition of two 



7 



Lorentz functions with variable weight, position and width as the input spectral function. 
The results are shown in Fig. 1 and 2a. Altough we are able to fit the numerical data with an 
accuracy better than 5 • 10~ 3 , several different fits are possible (only three of the innumerable 
possibilities are shown). Despite the chance of finding a spectral function fitting the data 
even more accurately, we do not think that it makes sense to decide between different fits, 
all having this extreme high accuracy (the G(t) data themselves are not exact due to the 
finite N and the finite resolution of the computer)! Therefore we conclude that structures 
with small weight like an additional peak at the Fermi level cannot be derived from the 
knowlegde of G{r). 

In Fig. 2b the result of a NCA calculation for the same set of parameters is shown. If the 
transformation (34) is applied to the spectral function the dotted line in Fig. 1 is obtained. 
We observe a significant difference between the NCA result and our numerical calculation 
(note the position of the maximum in the spectral function). Here we cannot decide which 
theory gives the better results because both are in a certain sense approximative (due to 
the finite N in the functional integral calculation and the leaving out of crossing diagrams 
respectively). 

Fig. 3a shows results for G(t) in case the f-level is equal to the chemical potential for 
different values of the hybridization. For V = the empty and singly occupied states have 
equal probabilities: 

n = n T = n i = - (35) 

This is the mixed-valence regime of the SIAM. With increasing hybridization more and 
more f-electrons are transferred to the conduction electron states. Therefore rif decreases 
with increasing V (Fig. 3b). This decrease is a little bit stronger in the NCA data but 
proportional to V 2 in both cases. 

We expect the spectral function to consist of only one Lorentz peak. With this input knowl- 
edge the fit procedure always leads to a well defined result (if the true spectral function really 
resembles a single Lorentzian) . Again the accuracy is very high but decreases with increasing 
hybridization. The results are shown in Fig. 4a. For zero hybridization we would get a delta 
function at uj = 0. If we increase the hybridization we observe a broadening and a shift of 
the peak to higher frequencies. Both effects are proportional to V 2 — this corresponds to 
calculations within a second order perturbation theory. 

The NCA spectral functions (Fig. 4b) give qualitatively the same picture. They are broader 
and shifted to higher frequencies than the functional integral results. 

The calculation of the f-susceptibility and the corresponding effective magnetic moment is 
straightforward within the theoretical scheme described above. We just add to the Hamilto- 
nian a term 

H B = -g^B (/|/ T - ///;) (36) 
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and calculate the f-occupancies raj and n^. The susceptibility is then defined by 



X = 



dB 



(37) 

B=0 



In Fig. 5 we show the dependence of the effective magnetic moment xl flg^B on the temper- 
ature for different values of the hybridization. For zero hybridization, the effective magnetic 
moment takes its maximum value at T = 0. For V > 0.1, due to the screening of the con- 
duction electrons, the effective magnetic moment decreases with decreasing temperature. Of 
course, the exact behaviour for T — > cannot be extrapolated from the data which are 
restricted to temperatures higher than k^T 0.2 . 



6 Conclusion 

In this paper have we investigated the U — oo single impurity Anderson model using a new 
functional integral technique in which the constraint of no double occupancy of the impurity 
site is fulfilled exactly by auxiliary complex fields. We obtained an effective action for the 
impurity as the starting point for the numerical calculation of the f-electron Green's function. 

We discussed the difficulties to extract detailed information for the spectral function out of 
the imaginary time Green's function. In the mixed-valence regime, where we expect that a 
single peak function is a good approximation for the spectral function, we can confirm a 
shift and broadening of the peak at the impurity level proportional to V 2 . 

Comparison with NCA results showed small but significant differences. The origin of these 
differences is not yet clear. 

Calculation of the susceptibility at the impurity showed the expected temperature depen- 
dence. While for zero hybridization, the effective magnetic moment reaches its maximum 
value for T — > 0, for finite hybridization we observe a decrease of xl Pd^B with decreasing 
T. 
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Figure captions 



Fig. 1 Green's function G(r) on the imaginary time axis (parameters: V = 0.22, £f = —0.2, j3 = 
3.0, D = 3.5, p — 1, N — 20 ). The crosses are the numerical result. Also shown are fits 
(nearly indistinguishable) calculated from the three spectral functions in Fig. 2a. The 
dotted line is the NCA result calculated from Fig. 2b. 

Fig. 2 a) Three different spectral functions whose corresponding G(t) (see Fig. 1) fit the 
numerical data equally well, b) Comparison of the NCA spectral function with fit 1 of 
Fig. 2a. 

Fig. 3a Imaginary time Green's function for Ef — p = 0, /3 — 2.0, D = 3.5, p — 1, N — 20 and 
different values of the hybridization. 

Fig. 3b Dependence of the f-occupancy n f on the hybridization (same parameters as in Fig. 3a). 

Fig. 4 a) Fits for the spectral function to the data shown in Fig. 3a with a single Lorentzian. 
b) The NCA spectral functions for the parameter set of Fig. 3a. 

Fig. 5 Dependence of the effective magnetic moment on the temperature for different val- 
ues of the hybridization (parameters: e { = —0.2,D = 3.5, p — 1,7V = 16). The line 
corresponds to the analytical result for V = 0. 
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